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rn ■ ABSTRACT 
0^ ' Recently, Madsen has argued that relativistic decays of massive neutrinos into lighter fermions and 



bosons may lead, via thermalization, to the formation of a Bose condensate. If correct, this could generate 
mixed hot and cold dark matter, with important consequences for structure formation. From a detailed study 
OO I of such decays, we arrive at substantially different conclusions; for a wide range of masses and decay times, 

we find that stimulated emission of bosons dominates the decay. This phenomenon can best be described as 
^ I a neutrino laser, pumped by the QCD phase transition. We discuss the implications for structure formation 

^ ! 

■ and the dark-matter problem. 

(N 

' Recent studies [1,2] suggest that a mixed dark-matter (MDM) universe could more readily account for 

, the large-scale power observed in galactic surveys [3] and recent COBE measurements of the microwave 

' background fluctuations [4]. These developments have stimulated attempts to provide a natural physical 

' , mechanism to generate MDM. We focus here on the novel and interesting idea of Madsen [5] that if a 



relativistic "heavy" neutrino, vh, decays into a light fermion, F, and boson, B, then a significant fraction of 
the bosons may form a condensate. 



^ ' Consider the processes: 



z/ff^F + B ;Ph^F + B , (1) 



where mp and ms are small enough (<J 100 eV) so as not to overdose the universe and mp^ms ^ rnn- 
Madsen argued that if the decay processes become effective at T = Td (when the time-dilated free decay 
time for a i^h of typical momentum is equal to the age of the universe) and if the vh are still relativistic, 
then the vh, F and B populations would come into thermal (but not chemical) equilibrium. Solving for the 
parameters of the thermal distributions and assuming that the abundances of B and F were zero at T^, he 
finds that, for a wide range of initial conditions, there is a Bose condensate. 

While highly suggestive, this heuristic argument is somewhat unsatisfactory. First, the process is not 
complete. In coming to the supposed equilibrium state only about half of the i'h^s decay. What happens to 
the remaining vh's is unclear. Second, a sudden equilibration at may be an inadequate description. To 



address these issues we have solved the Boltzmann equations for the reactions (1). 

We find that if the decay rate is sufiiciently high (T^ > sJmH'mB) then the equihbration will be preceded 
by a burst of stimulated decays into very cold (p <C T) bosons. This is an exponentially unstable process 
i.e. a neutrino-laser. As a result of this non-thermal process, close to half of the vr'^ decay. The resulting 
vh and F distributions have a 'grey-body' form and the cold B's are also non-thermal, with an abundance 
well in excess of the thermal-equilibrium Bose-condensate. 

Allowing only the reactions (1), the evolution of the occupation number distributions /j, i =H,F,B, are 
described by the following Boltzmann equations: 

h = {himlTo/moEiPi) J dEj^i[fH{l - + /s) - /s/f(1 - /h)] (2) 

where fi is the time derivative of fi at a fixed comoving momentum, /it = — 1,/if = 1 and hs = 2(1) if 
B = (^)B, Pi is the physical three-momentum, and mo/2 is the three-momentum of the decay products in 
the vh rest-frame with tuq = m\ — 2{m\ + m%) + {m\ — m?pY /m\. Fq is the free decay rate for a vh at 
rest. The integration is over the energy- conserving plane Eh = Ep + Eb, with limits on Eb,Ep 

E^ = {mo/2mH){EH^l + ^ml/ml ± ph) (3). 

(fc = B,F) as shown in figure 1. These limits follow from purely kinematic considerations. For high Eh, F 
and B both come off in the forward direction with similar energies. For 

Eh = E^ = {m% + m% - m%)/2mB — m%/2mB 

zero-momentum bosons are accessible. If the final states were empty the uh's would decay into products 
with energies distributed uniformly on the range (3). Equation (2) then follows straightforwardly from the 
inclusion of the quantum mechanical statistical weights for the forward and inverse decay processes. 

To set the initial conditions we assume that at some early time vh, F and B were all in chemical and 
thermal equilibrium, but that the F and B decoupled prior to the QCD phase transition at T w 100 MeV 
while the vh did not; this specifically requires F fe,t'ju. Following the phase transition the temperature 
of the uh's was increased by a factor rj^/^ ~ 2 relative to the F's and B's, where r] is the ratio of statistical 
weights before and after the phase transition. Equations (2) can be integrated numerically, but we can gain 
insight into the general properties of the solution by considering some limiting regimes. 

Consider the equation for fn- Since 1 — /f < 1 and the range of the integral is E^ — = mopH/mn, 
the rate for terms independent of Jb is ~ mH^o/En- This is just the time dilated free decay rate for a 



typical energy uh and therefore will be negligible for T = min(To, (Tomij)^/^), where Tq is such that 
-ff(To) = Tq. The terms involving fs, however, can be important long before T^. Neglecting the other terms 
we have 

fB = {hBTomlfB/moEBPB) J dE{fH - fp) (4) 

so fs will grow exponentially, driven by any initial imbalance between the vh's and F's. The limits on the 
integral in (4) are E^(Eb), which are given implicitly by Eq (3). The prefactor 1/EbPb suggests the growth 
rate will be greatest for the lowest momentum bosons which are accessible. Consider first the case T 3> -E* 
which must hold at sufficiently early times. The minimum boson energy for a typical energy uh {Eh ~ T) 
is E^ {mB/TnH)^EH. Setting Eb,Pb ^ E^ and since the integral in equation (4) will be on the order 
Eh ~ T we obtain 

fs ^ {Tom%/m%T)fB (5) 

so the growth rate is Fiase — ro(m_f//m_B)'^(-E*/T). This increases with time while the expansion rate 
decreases, so these will be equal at 

Tiase ^ (m^T2/m|)i/3 (6) 

We have assumed here that T ^ E^. For T <^ E^ we find that the growth rate decreases faster than the 
Hubble rate so if the process is not effective at T ~ then it never will be. The condition Tiase ^ E^ sets 
a lower limit on the decay rate, or equivalently on the decay temperature: T4 > ^JruHrfiB. 

What we have here is a neutrino laser; the low momentum boson occupation numbers will grow expo- 
nentially via stimulated decays, feeding off any initial imbalance between the z/h's and F's, and terminating 
when this is driven to zero. The momentum of the bosons produced is pB — ('Tis/"^!/)^) so for a large mass 
ratio these will be very cold compared to the typical thermal energy. 

The lasing will result in fn = fp = {fn + fp)/'^ (superscript indicates initial values) which is a fermion 
analogue of a 'grey body' spectrum, and the total number of cold bosons produced is just (uh" — ^^F°)/2 = 
(1 — T]~^)nH° /2. The number of hot bosons is just the initial number ub" = 4n^/(3ry/iB). This then gives 
the fraction of cold bosons after lasing to be ~ 80(65)% for B = (^)B. The result of lasing is therefore 
qualitatively similar to the thermalization calculation of Madsen, but with important differences: The lasing 
process occurs much earlier than thermalization; the initial hot bosons are unaffected and the decays occur 
exclusively into the cold component. In the equilibration calculation, only about half of the decays go into 
the cold phase, and the cold fraction is ~ 42(26)%, about half the yield from lasing (these numbers are 
slightly different from those calculated by Madsen as we have allowed for the finite initial boson abundance). 



After lasing the temperature will eventually fall to T^. What happens then depends T4. For < 
the low momentum bosons produced by lasing will be effectively decoupled (they can only be reached from 
Eh El, ^ T and so this coupling will be exponentially suppressed). The hot bosons and fermions will 
equilibrate, but they do this without generating any further bosons. For T4 > E^, on the other hand, the 
laser generated bosons are still accessible from typical energy decays, so we expect that some of the cold 
bosons will be reabsorbed and that the post-T^ distributions will adopt the form predicted by Madscn. 

In either case the decay process is not yet finished since there is still a finite vu abundance. The 
remaining vrs (58% of the original abundance) decay into hot B's and F's once T ^ min(mii-, T^). For 
lasing, the final fraction of cold bosons is {l — 'r]~^)/{2 + 8/3r]hB) — 37(34)%. For equilibration (T^ > Et,) the 
numbers are 19(13)%. When and how the final decays occur depends on T^. For T^, > Et,, the populations 
equilibrate at and are then stable until T ~ niH when the heavy neutrinos go non-relativistic at which 
time they decay. For mjj <Td < Et, we expect the hot particles to come to thermal equilibrium at T ^ 
without generating any further cold bosons. The remaining i/jj's decay at T ~ tuh, with no subsequent 
equilibration. For y^mnmB <Td < ttih no equilibration occurs and the particles decay at Td when they are 
non-relativistic. The last possibility is < ^Jm^jriiB in which case no cold bosons are produced and the 
vh^s decay into very hot products at T^. 

These predictions are confirmed by our numerical solution of (2). In Figure 2 the evolution of the 
occupation numbers are plotted for masses and decay temperature chosen to illustrate the lasing phenomenon. 

In all of this one must be wary that the new fermionic degree of freedom F does not take up its 
full statistical weight prior to the weak interaction freeze-out, thereby violating the bounds from standard 
big bang nucleosynthesis [6]. (The B's are not of concern if the stimulated decay is dominant since they 
are very low energy.) We therefore require Tiase < 2.3 MeV since the number changing reactions for the 
neutrinos are decoupled below this temperature. From equation (6) and the definition of Td this requires 
Td < (2.3MeV)3/2j„|^-5/2 j.^ ^ j.^ ^ 2.3MeY{mB/mH)^^^ otherwise. The former, together 

with the lower bound Td > ^JmuraB sets an upper limit niH ^ ^/2MeVmB ~ lOkeV for a boson mass which 
would close the universe. These, together with the requirement that mn > itibjITIf, are shown in Fig. 3, 
which displays the region of Td — niH space where cold bosons are produced. 

The fraction of cold particles is important for structure formation, and provides a firm (modulo factors 
such ashB,T]) prediction of the lasing model. Also important, however, is ms and the momentum distribution 
of the hot bosons, as they will determine the maximum Jeans mass [7]. The attraction of the 70/30 cold/hot 
ratio MDM model [1] is mostly that the low hot-fraction reduces the mass of the neutrino by a factor ~3, 
therefore increasing the maximum Jeans mass by a factor of ~10. Here the final abundance of bosons is just 
{l + 4:/3r]hB)n'^. If the bosons alone constitute the dark matter then their mass will be only slightly less than 



the standard value. However, the boson momentum distribution will also differ from the standard value. For 
Td > -B* the mean momentum of the hot bosons is almost identical to the standard value so the maximum 
Jeans mass will be very similar to that in the standard HDM model. This may be problematic. In the case 
of lasing with Tj, > mjj the hot bosons are ~ 40% cooler than in standard HDM making matters even worse. 
However, for the lower triangle ^mnmB <Td < m^, the temperature of the hot bosons can be larger than 
the ambient temperature by a factor ^ niH /Td- This can be as large as ^ y/ran /mB ^ (2MeV/mB)^/^ ~ 20. 
It is also possible that the fermion and boson masses are similar, or that other particles, e.g. the V/j, decay 
into the same products; in either case this would increase the maximum Jeans mass. 

We note that a number of variations on the preceding discussion exist, which may have different quan- 
titative predictions. For example, if vr = ^n, and is heavier than Ut, then the QCD phase transition 
would pump the v^j^ — vh system, because charged-current production of Ur is energetically disfavored; or 
the uh, F and B could all be new particles outside the standard model. Investigation of such possibilities is 
beyond the scope of the present work. 

We have thus far only considered the scenario where Fy^ Ve^^n- Consider briefly the scenario where 
the light fermion is one of the known neutrinos, and is therefore coupled at Tqcd- Madsen's calculations 
[5] indicate that in this scenario too, Bose condensation may occur. With regard to our calculations note 
that there is no initial fermion imbalance to drive the lasing. In fact, we conclude that direct solution of the 
Boltzmann equations (Eq. 2) results in no condensation. 

In conclusion: we have shown that for decay temperatures Td > s/mHtriB that thermalization is in- 
evitably preceded by a runaway stimulated decay into cold {p ~ {mB/'mH)^T) bosons. This 'neutrino lasing' 
provides a relatively natural way to generate MDM. For decay temperatures Td » m%/mB the populations 
will thermally equilibrate, some of the cold bosons will be reabsorbed, and the final cold boson fraction is 
~ 15%. For Td <C m%/mB the laser generated bosons survive to the present and the cold fraction is ~ 35%. 
If the boson alone constitutes the dark matter, then thermalization results in a maximum Jeans mass very 
similar to that in the standard HDM model; in contrast, the lasing scenario allows the interesting possibility 
of a higher Jeans mass. 

We greatly appreciate numerous discussions with Arif Babul, Dick Bond, Malcolm Butler, Savas Di- 
mopoulos. Lev Kofman, Lenny Susskind, Chris Thomson, Scott Tremaine and Ned Wright. 
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FIGURE CAPTIONS 

1. Allowed Eb for decay of uh with energy Eh (logarithmic scales). Of particular importance for our 
analysis is the form of E^ since the lasing occurs into the lowest momentum states which are accessible. 
Low energy vh's decay into high energy products. Particles with Eh = -B* 7X1% /niB can decay into zero 
momentum bosons and for Eh > -E*, Ej^ ~ {mB/mH}^EH- 

2. Upper panels show p^f, the number of particles per log(p) for three output times: the initial time, after 
lasing and after all vh's have decayed. The masses were tuh; mp; thb = 1; 0.03; 0.03 and = 1. The 
lower panel shows the total number of bosons as a function of time. The lasing epoch, which injects cold 
bosons exclusively, followed by a trans-relativistic decay epoch, injecting hot bosons, are evident. Most 
striking is the bimodality of the boson distribution which is clearly divided into hot and cold components. 

3. Allowed regions for generating mixed dark matter (logarithmic scales). In the upper triangle lasing 
occurs, but is then followed by equilibration resulting in a reduction in the number of cold bosons. In the 

lower triangles lasing occurs. Above the dashed line, the sequence of events is lasing; equilibration of the 
hot components; and finally decays at T ~ mn- Below the dashed line, the particles become non-relativistic 
before the decay process terminates. 



